N76-15029 


AIPSHIP STRESSES DUE TO 
VERTICAL VELOCITY GRADIENTS 
AND ATMOSPHERIC TURBULENCE 

Duncan Sheldon* 


ABSTRACT ^ Monk's potential flow method is used to calcu- 
late tho resultant moment experienced by an ellipsoidal 
airship. This method is first used to calculate the moment 
arising from basic maneuvers coisidered by early designers , 
and then expended to calculate the raoment arising from ver- 
tical velocity gradients and atmospheric turbulence. This 
resultant moment must be neutrali zed by the transverse 
force of the fins. The results ihow that vertical veloc- 
ity gradients at a height of 6000 feet in thunderstorms 
produce a resultant moment approximately three to four 
times greater than the moment produced in still air by 
Idealistic values of pitch angle or steady turning. Real- 
istic values of atmospheric turbulence pioduce a moment 
which is significantly less than the moment produced by 
maneuvers iu still air. 


INTRODUCTION 

At one time airship design was a highly organized and systematic 
activity, and hundreds of papers have been writt'^n on the subject. 

The period of greatest activity was from 1910 to .938. However, in 
spite of careful efforts several notable disasters occurred. Some 
were at least partly the result of political considerations; examples 
are the American ship Shenandoah and the British 3hip R- 1 0 1 . The 
most spectacular of all, the Hindenl urg disaster, was of course due 
to the use of hydrogen as lifting gas. With the exception of the 
inadvisable use of hydrogen and the deterioration of the hull of the 
R ^lQl » most well-].nown dirigible disasters were connected either 
with atmospheric turbulence or vertical wind currents in storms or 
above mountains^. 
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The British ship R-38 buckled in the middle and broke in two because 
of a strong wind gust (1922) . At the time the airship was already 
experiencing significant stresses arising from a sharp turning 
maneuver. The Shenandoah perished in a 70 mph squall (1924). As the 
result of a navigational error the Akron was drawn into the center of 
a storm. While maneuvering upward to offset a downdraft, its lower 
rudder hit the ocean and the airship fell into the sea (1933). The 
Macon lost its top rudder during a squall and was also lost at 
sea (1935). 

As a result of the R-“38 crash the Royal Aeronautical Society estab- 
lished the R-38 Memorial Prize. In response to this competition 
three exceptionally detailed airship design papers were published^* ^ . 
This was in 1923, and taken together they constitute probably the most 
detailed airship design analyses available in English. Most later 
work was a refinement of methods discussed in these articles. One can 
even view the design of the Graf Zeppelin and her sister ships (1928- 
1938) within the context of the methods presented by these British and 
American authors. Of course the principal ingredient missing from 
these relatively early papers is the practical experience and 
full-scale data obtained by the German designers. However, there were 
no basic changes in the relevant technology in the years from 1923 to 
1938. 


An important part of the early design work was the highly ingenious 



Kelvin impulses. Under most conditions Munk's theory is in surpris- 
ingly close agreement with full-scale experiments^ - 


7-9 

As pointed out in several recent articles , the technology relevant 
to airship design has undergone an extraordinary expansion along with 
all other aerospace activity. Modern computers and modern knowledge 
of structural dynamics permit analyzing the airship's structure as a 
whole. It is essential to apply our current knowledge of atmospheric 
turbn.lence and vertical wind currents to these structural calcula- 
tions. Safety is the overwhelming design consideration applicable to 
future airships, and relating atmospheric hazards to structural 
integrity holds the greatert promise of assuring safe operation. It 
might be argued that damaging atmospheric effects can usually be 
avoided, particular!’- during non-scheduled flights. The record of the 
German pilots server to establish this to some degree. But the impor- 
tance of scheduled operations also requires that atmospheric hazards 
be given careful consideration. 


The purpose of this paper is to show how our present knowledge of 
the atmosphere can be combined with Munk’s equations to calculate 
the resultant moment or. an airship arising from vertical currents 
and atmospheric turbulence. Approximate results are given for the 
resultant moment experienced by a 1,000 foot long ellipsoidal air- 
ship with a fineness ( length- to-diameter) ratio of 5. This is the 
shape suggested for a ’’basic” airship considered in a recent design 
study by Mowforth**-^ . These results are compared with the moments 
arising fmm pitch angles and steady turning rates in still aii 
which were taken into consideration by the early designers. 
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AERODYNAMIC CONSIDERATIONS 


Munk's Equations 


The motion of airship hulls gives rise to an air flow that is well 
approximated by potential flow. There may be a large resultant moment 
of the aerodynamic forces r but only a comparatively small lift and 
drag. With wings the conditions are different as there is consider- 
able lift- Since the momentum of the flow is not necessarily in the 
direction of motion of the hull, a principal axis problem presents 
itself. Strictly speaking, we should distinguish between the momentum 
of the flow and the Kelvin impulse of the flow, but Munk himself dis- 
regarded this difference and we have no need to make the distinction 
here. The net resultant moment is expressed in terms of the 
volumes of the apparent additional masses of the hull. The apparent 
additional mass of a solid moving through a fluid along one of its 
principal axes is simply a proportionality constant expressing the 
resistance to accelerations along the axis offered by the fluid 
itself. Note that it is not a measure of the inertia of the solid, 
because the solid need not have any mass at all. In this case all of 
the energy is stored in the flow, and the apparent additional masses 
along each principal axis are equal to the apparent masses. The 
effect of the fluid surrounding the solid is, however, fully described 
by assigning to the solid an apparent additional mass in addition to 
its original or actual mass . The apparent mass of a circular cylinder 
in a uniform two-dimensional stream is PTfr^; and for a sphere in a 
three-dimensional uniform stream its value is 7(^trr^)p. Here r and p 
are radius and density. Apparent volume is obtained from apparent 
additional mass by dividing by the density. 




' % 


Munk shows that an airship hull, flying steadily under an angle of 
attacked and with the velocity of flight V experiences a resultant 
couple of the magnitude 


M = (K2--K]^i sin 2 << 


( 1 ) 


where K, and K2 denote the apparent volumes with respect to the 
longitudinal and transverse principal axes of the hull. This moment 
is unstable, consequently fins are required for stabilization. Munk 
also calculates the transverse force on an airship (with circular 
cross section) turning under an angle of yaw: 


dP = dx[(k,-k, )|2 -^sin2(}, 

2 I dx 2 

+ k'V^ -oor^(t, + k'V^ 3 X — cos^4i] 

R K dx 


(2) 


where 


dp = Transverse force acting over a differential 
length along the longitudinal axis 

dx = Differential length along the longitudinal axis 

k]^ = (Hull volume) 
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k 2 = (Hull volume)/K 2 

k* = Ratio of the apparent hull moment of inertia 
about the aerodynamic center to the moment of 
inertia of the displaced air 

X = Position on the longitudinal axis relative to 
aerodynamic center 

S = Area of circular cross section at x 
(J) = Yaw angle 
R = Turning radius 
V = Airship velocity 
p = Density of air 

This expression of course does not contain the air forces on the fins. 
Munk's theory also yields a closed form expression for the pressure 
distribution over any ellipsoid inclined at an arbitrary angle to the 
flow- The first term on the right-hand side of Equation 2 can be used 
to calculate the longitudinal distribution of forces resulting from a 
vertical gust. In this case the yaw angle in Equation 2 is identified 
with the angle of attack 




4> = t.'n ^ (3) 

where u is the transverse velocity and V is the forward velocity. 

Munk assumes the airship has a variable effective angle of attack 
along its axis. The magnitude of the superposed angle is tan~^ (u/V) , 
where u generally is variable. The momentum produ ed at each portion 
of the airsiiip is the same as the air force at that portion if the 
entire airship had that particular angle of attack. Consequently, 
Equation 2 can be used to determine the moment experienced by an air- 
ship as it moves through a vertical velocity gradient. In this case 
we assume the pilot is able to hold the airship on a straight course 
in inertial space without yaw or pitch. Equation 2 will also be used 
to calculate the moment resulting from a turning maneuver. Equation 1 
provides a direct method of calculating the bending moment when the 
only disturbing force is due to pitch. 

Moment Response Function 

Munk*s theory can be extended to calculate the transverse forces caused 
by atmospheric turbulence. It is assumed the pilot is able to hold the 
airship on a straight course in inertial space without yaw or pitch. 

We begin by attributing to a circular cross-section of area S the 
virtual mass Spdx just as if the cross-section were part of a circular 
cylinder immersed in two-dimensional flow. The transverse force acting 
on this cross section as a result of the velocity perturbatio \ u - 
UqC^^^ is 

f = j>Sdx(io)uoe^'^^ (4) 
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*o = 2irf = 2ir^ = kpV (5) 

where a) is the. angular frequency of the perturbation, f is the cyclical 
frequency, V is the forward velocity of the airsh , X is the wave- 
length in the forward direction, and kp is the propagation constant for 
a particular wavelength. It is convenienc to take the geometric center 
of the ellipsoid as the origin of our coordinate system. Then the 
moment experienced by the airship, per unit velocity perturbation, is 
given approximately by 




+L/2 

\ S(x)e^ 

^-L/2 


(u)t+kx) 


X dx 


( 6 ) 


where L is the length of the airship. 

Uniform S is not a candidate hull shape, but this case leads to the 
simplest form of the moment response function. If S is uniform, 
the result is 




SL ( 


sin(irC ) 


cos (^5) ) 


(7) 


where 5 = l/A . This is the long wavelength approximation, and 
approaches zero as $ approaches zero. For short wavelengths, 
* 3 ^ 1 , the bending moment at the longitudinal positions of maximum 
transverse velocity is the important consideration. In this case 

(« I . i£V J - 1, (8) 

I Uo‘ kp ^ 

For an ellipsoidal airship with a fineness ratio of 5 we set 

S(x) = ^ |(L/2)2 - 


(9) 


and use Equation 6 to obtain 

TT ^3 


1 u'ol' 50 



A 

3 cos k 




A 

3 sin Jc 




(10) 


where k = Lk /2, again this is the long wavelength approximation, and 
the right-hand side of Equation 10 tends to zero as ^ tends to zero. 
The short wavelength approximation, Equation 8, still applies 
provided 


ds A 

5x S 


« 1 


( 11 ) 


161 


y 



I 




) ‘ ^ 
■4 

i 

h 




METEOROLOGICAL CONSIDERATIONS 
y t i nd Gradients in Thunderstorms 

Using Munk * s equations we can calcuJate the force at each section of 
the hull of a representative airship for vertical wind currents known 
to exist in the atmosphere. Typical values for a thunderstorm are 
considered. Taken together Figures 1-3 enable us to obtain a good 
approximation of the vertical currents and horizontal scale of thunder- 
storms. Figure 1 shows information which was obtained to describe tne 
atmospheric effects an airplane experiences as it flies through a 
thunderstorm. Figure 2 suggests that the vertical velocity profile 
given in Figure 1 is applicable above about twenty thousand feet. 

With the help of Figure 3, we can construct a similar thunderstorm 
profile for an altitude of approximately 6,000 feet, which is a 
typical operational altitude. These diagrams give no information 
about the severity of the turbulence; they can only be used to study 
the airship stresses arising from vertical currents. However, we can 
observe that the region of "violent turbulence" extends much further 
than the region of severe vertical currents. Figure 1 represents the 
vertical velocity profile in the plane of travel of the thunderstorm; 
the updraft usually has a fairly uniform cross section of about 10 
miles traverse to its line of travel. Figure 3 shows that below about 
ten thousand feet the vertical flow is not quite as constricted as at 
higher altitudes. Let us therefore make the approximation that at 
6,000 feet the horizontal scale of the currents is twice as large as 
shown in Figure 1. Accordingly, the magnitude of the vertical current 
distribution is cut in half (due to the transverse extent of the storm 
the flow is treated as two-dimensional). This means that the more 
severe vertical velocity gradients in the horizontal direction are of 
the order of 0.2 ft. /sec. /ft. 

Atmospheric Turbulence 

It is possible to describe the energy distribution of atmospheric 
turbulence as a function of wavelength by the same relatively simple 
formula for the following three forms of turbulence 

(i) clear air turbulence near the ground, 

(ii) turbulence near and in cumulus clouds, 

(iii) thunderstorms . 


The distribution of energy density S^(k) at diffeient wavelengths 
(1/k) of the vertical component of the turbulence moy be given by 


where 




1-1- (8/3) (Lik)^ 
jl + (Ljk)2 J 


= Root mean square vertical velocity 
Lg * Scale of turbulence 


( 12 ) 


= 1.339 (2eL8) 

for each patch of turbulence. Twice the total energy per unit mass of 
air equals the mean square of the turbulence velocity so that 
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(13) 


Thus/ twice the total energy per unit mass of air is given by the area 
under the curve of kS^(k) against log^ k and the area under the curve 
kS^(k)/o^ is unity. Usually Equation 12 is adjusted to fit experi- 
mental data by selecting Lg so that the calculated and experimental 
distributions kS^k) have a maximum at the same value of k. A compari- 
son of theoretical and experimental distributions is shown in 
Figure 4 . We shall follow the common practice of referring to 

S^(k) as a power spectral density even though in reality it is a mean- 
square- value-density spectrum, 

NUMERICAL RESULTS 

Wind Velocity Gradients 

The resultant moment experienced by the airship is evaluated from 

,+L/2 

AT? 

(14) 




L/2 


(g)xdx 


by using Equations 2 and 3 and setting u = (du/dx) {L/2 - x) . The 
results are shown in Table I for (da/dx) = 0 . 01 , 0.1, 0.2, and 0.3. 

Atmospheric Turbulence 

If H(k) is the response function describing the resultant couple when 
the airship is subjected to a transverse velocity wave of unit 
amplitude, then the mean square value of this moment is given by 


V / 


n,+0O 




r.m.s. 


|H(k)\ S^,(k)dk 


(15) 


A A 

if S^(k) is a stationary function. Using dk ~ dk/(1TL) and dk = 
kddoggic) Equation 15 becomes 


M 


r.m.s. TTL 


After setting 


lH(k)|^ S^(k)1cd(log 


(16) 


'-00 


H(k) 


M 

Ue 


(17) 


Equation 16 was used to evaluate in response to the atmospheric 

power rpectral density function given'by Equation 12. Two cases were 
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considered : (1) Captive balloon data taken at a height of 1,000 

feet with -• 3-54 fps, and (2) Data obtained from an airplane flight 
at 40,000 feet in a thunderstorm. The results are shown in Table I. 


Summary 


The resultant moments obtained by the various methods discussed in 
this paper are compared in Table I. Dashes are used where an entry 
is not applicable. These moments are an important measure of the 
airship's stress because they must be neuvralized by the transverse 
force of the fins. The first five cases, which include rectilinear 
motion at a constant pitch angle and steady turning without pitch, 
are conditions in still air which were considered by the early 
designers. The angle of yaw corresponding to the turning radius R 
was obtained by Munk^ 


4 


. L 1 

2R k2 - k]^ 


(18) 


Equation 1 was used in Case 1, and Equation 2 was used in Case 2. 
Agreement of these two cases serves as a check on the numerical 
methods and also confirms, with remarkable accuracy, the approxima- 
tions Munk used in deriving Equation 2. 

Cases Six through Eleven correspond to situations where our current 
knowledge of the atmosphere was used. When a uniform vertical veloc- 
ity gradient was considered, the vertical velocity was assumed to be 
zero at the tail and increase in the direction of flight. The resul- 
tant moment for Case Nine is less than Case Eight because the sine of 
the angle 2^> contained in Equation 2 decreases as 4> increases beyond 
45®. The data for the thunderstorm were obtained at 40,000 feet and 
are not fully satisfactory for our purpose. However ^ the density was 
adjusted to this height, and the result corresponding to a direct 
application of the power spectral density equations is included. These 
results show that the values of atmospheric turbulence found in the 
literature produce a moment which is significantly less than the 
moment produced by realistic maneuvers in still air. However, the 
vertical velocity gradients et an altitude of 6,000 feet in a thunder- 
storm produce a moment which is three to four times larger than the 
moment produced by maneuvers in still air. 
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Coefficients of Additional Mass of Ellipsoid: k^ = 0.059, k 2 = 0.895, k* = 0.701 
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